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Preface 


This  report  is  the  result  of  questions  raised  about  variations  in  the  recorded 
values  of  temperature  of  air  passed  through  balloon-borne  atmospheric  sampling 
devices.  The  results  of  this  investigation  provide  essential  information  to  the 
scientific  users  of  the  system  under  examination,  and  have  implications  *o  users 
of  balloon-borne  systems  where  external  temperatures  are  a critical  parameter. 

The  idea  for  this  report  was  generated  in  conversations  with  Mr.  Thomas 
Ashenfelter  of  the  National  Oceanic  and  Atmospheric  Administration-Air 
Resources  Laboratory,  Mr.  John  Ground  of  New  Mexico  State  University,  and 
Captain  J.  Robert  Greenlee  of  the  Air  Force  Geophysics  Laboratory  (AFGL). 

The  author  is  indebted  to  his  colleagues  at  AFGL  for  support  and  guidance  in 
preparing  this  document:  to  Mr.  George  Nolan  for  suggesting  the  present  format, 
to  Mr.  James  Dwyer  and  SSgt  Robert  Dumont,  respectively,  for  checking  the 
calculus  and  computations,  and  to  Mr.  Ralph  Cowie  and  Mrs.  Catherine  Rice  for 
reading  the  initial  draft  and  providing  many  thoughtful  suggestions.  Special  thanks 
go  to  Mr.  John  Ground  for  rapidly  forwarding  the  data  in  Tables  11-13,  and 
Captain  J.  R.  Greenlee  for  researching  flight  records  and  processing  the  large 
volume  of  data  which  is  so  succinctly  presented  in  Table  14. 
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An  Examination  of  the  Temperature  Measuring 
Accuracy  of  a Flowmeter  System  Used  With 
Balloon-Borne  Atmospheric  Samplers 

I.  INTRODI  CTION 

Every  year  numerous  balloon  missions  are  conducted  to  carry  samplers  to 
predetermined  altitudes,  usually  between  15.  2 km  (60  k ft)  and  36.6  k m (120  k ft), 
to  gather  specimens  of  atmospheric  constituents  for  various  scientific  groups. 
Motor  driven  fans  or  nitrogen  driven  aspirators  draw  air  through  various  types  of 
filters  which  remove  gaseous  or  particulate  constituents.  One  sensor  monitors 
the  flow  rate  and  temperature  at  the  inlet  and  another  the  exhaust  of  the  sampler. 
This  sensor  is  a 12.7  cm  (5  in.)  diameter  cylinder  about  15.2  cm  (6  in.)  in  length. 
Temperature  and  flow  rate  are  monitored  independently  by  a thermistor  and  turbine 
driven  chopper  respectively.  These  two  measurements  are  used  to  calculate  the 
volume  of  air  passed  through  the  filter  element.  The  accuracy  of  the  volume 
calculation  is  directly  proportional  to  the  accuracy  of  the  data. 

Over  the  last  several  years  interest  has  been  growing  in  two  characteristics 
of  the  temperature  data.  First,  during  a typical  flight  the  measured  temperature 
values  tend  to  oscillate  about  some  reference  temperature.  Second,  a large  differ- 
ence in  average  temperature  between  flights  at  similar  altitudes  has  been  observed. 
The  purpose  of  this  analysis  is  to  see  if  these  results  are  being  caused  by  the 
instrumentation  or  if  the  instrumentation  is  correctly  reporting  the  situation. 
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Each  of  these  characteristics  is  either  real  or  apparent.  Real  is  defined  as 
an  actual  change  in  air  temperature;  apparent  is  defined  as  no  actual  air  tempera- 
ture change.  Change  is  defined  as  variations  of  the  temperature  indicated  by  one 
system,  or  the  difference  in  average  temperature  between  two  systems.  This 
analysis  will  explore  the  mechanisms  which  can  cause  apparent  temperature 
changes. 

Causes  of  apparent  temperature  changes  are  either  internal  or  external  to  the 
system.  Internal  causes  are  the  manufacturing  tolerances  of  the  system  compo- 
nents. These  tolerances  result  in  the  overall  system  accuracy  which  corresponds 
to  the  error  envelope  of  the  data,  and  the  precision  of  an  individual  piece  of  equip- 
ment. The  only  external  cause  for  apparent  temperature  change  which  has  been 
identified,  is  radio  frequency  energy  from  the  onboard  data  transmitter.  Labora- 
tory testing  has  verified  that  as  the  thermistor  resistance  increases  to  the  1 meg- 
ohm region,  energy  radiated  from  the  transmitter  can  affect  the  temperature  data. 
The  original  thermistor  has  a resistance  characteristic  which  exceeds  1 megohm 
at  about  -60°  Celsius. 

Reflecting  on  the  above  facts  in  light  of  the  availability  of  curve -matched 
thermistors  having  lower  resis-  mce  and  smaller  mass,  the  author  decided  to  test 
a different  thermistor  while  cr  .ically  examining  the  original  thermistor  circuit. 
The  following  paragraphs  docu  nent  the  preliminary  calculation;:,  laboratory  tests, 
and  flight  data  on  which  the  conclusions  are  based. 

2.  TIIK  KLORMETER  SYSTEM 

V 

The  PR-3  Flowmeter  System,  * was  purchased  by  the  Department  of  the  Air 
Force  at  Goodfellow  AFB,  Texas  in  June  1969.  Temperature  is  sensed  by  a 
thermistor  located  in  the  air  flow  and  recorded  incrementally  as  a varying  direct 
current  from  zero  to  100  pA. 

In  this  report  the  term  ''orig!nal  system"  pertains  to  the  original  PR -3  circuit 
with  the  Fenwal  rod  thermistor  QB41J3;  the  term  "modified  system"  pertains  to 
the  author's  circuit  as  described  below  with  the  Fenwal  bead  thermistor  UUB31J1. 


1.  PR -3  Flowmeter  System  Technical  Handbook  (1969)  ASD  document  No.  3330 
by  the  Applied  Science  Division  of  Litton  Systems,  Inc.,  for  the  Dept,  of 
the  Air  Force  under  Contract  No.  F41616-69-R-0442. 
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2.1  The  Original  Circuit 


2.1.1  ACCURACY 

Figure  1 portrays  the  critical  parts  of  the  circuit  under  discussion.  Wire  and 
contact  resistances  are  small  enough  to  be  negligible.  From  Ohm's  law. 


t + r4  + r32 


where  the  subscript  t denotes  the  theoretical  current  which  is  dependent  upon  the 
values  of  ideal  resistors  and  an  ideal  source  E.  If  the  components  are  allowed  to 
vary  according  to  the  manufacturers  specifications,  this  same  expression  can  be 
used  to  define  the  real  current  I , where 


Tr  = Ra  + R4  + r32 


KEY 

D1 

1N4736A,  6.8V + 5% 

R4 

60.4k  ohms  + 1% 

Q1 

2N1308 

R32 

thermistor  QB41J3  10k  ohms  + 1%  at  25°C 

ra 

4600  ohms,  assumed  + 0%  see  K 

K 

ammeter  accuracy  + 2%  of  full  scale 

Figure  1.  Original  Temperature  Sensing  Circuit 


See  PR-3  Flowmeter  System  Technical  Handbook,  drawing  number  D-17788- 
236465,  page  37. 


To  include  the  effect  of  the  ammeter  which  displays  and  records  the  current, 
measured  current  is  expressed  as: 


t _ t 4.  

m r ^ " R.  + R„  + R 


+ K 


32 


(3) 


where  K is  the  full-scale  error  of  the  meter  and  recorder  (Stout  ). 


The  absolute  error  (U.  S.  Department  of  Commerce  ) of  I is  AI  = I - L. 

r m m m t 


Taking  Im  as  the  dependent  variable,  AIm  is  given  by: 


81 


81 


81 


81 


81 


AIm  ' 8K  AK  + 8E  AE  + 8R, 


AR.  + 


AR . +■ 


AR, 


A 8R4  4 8R32  32 


(4) 


This  analysis  initially  will  express  percent  error  of  current  over  a tempera- 
ture range.  Terms  E,  R^,  and  R^  are  independent  variables.  Resistor  R^^  is  a 
function  of  temperature,  and  K is  the  meter  accuracy  which  is  defined  by  the 
manufacturer  as  ±2  percent  of  full  scale.  The  effect  of  K will  be  considered  inde- 
pendent of  R32  and  related  to  current  later.  Term  K includes  all  the  meter 
inaccuracies;  therefore. 


ara  = ° 


and 


3K 


8Ra 


= 0 


(5) 


The  error  of  Rg2  is  a function  of  temperature  and  will  be  introduced  when 
numerical  values  are  calculated.  Equation  (3)  will  be  differentiated  to  obtain  the 
terms  necessary  for  Eq.  (4). 


81 


8K 


8 

E 

+ K 

(6) 

‘ 8K 

ra 

+ R4  + R32 

8 

E 

+ 9K 

(7) 

8K 

ra 

+ R4  + R32 

+ 5TT 

= 0+1 

(8) 

= 1 . 

(9) 

2. 


Stout,  M.  B.  (1950)  Electrical  Measurements,  Prentice  Hall,  Inc.  Pages  405 
to  408  have  a good  discussion  on  errors  in  electrical  indicating  instruments. 

3.  U.  S.  Department  of  Commerce  (1967)  Handbook  of  Mathematical  Functions. 

NBS/Applied  Mathematics  Series  No.  55,  page  14,  for  a discussion  of  error 
analysis. 
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And  for  the  regulated  voltage,  E 


91  a 
m _ 9 

9E  9E 


RA  + r4  + K32 


+ K 


9 E ^ 9K 

9E  ' RA  + R4  + R32  TO- 


Ra  + r4  + r32 


+ 0 


RA  + R4  + R32 


For  the  first  of  the  resistors. 


91 


m 9 


9Ra  9Ra 


RA  + R4  + R32 


+ K 


9K 


9ra  Ra  - r4  + r32  3RA 


0 - E • 1 


<RA  + R4  + R32* 


+ 0 


-E 


(RA  + R4  + R32) 


2 * 


Similarly  for  the  other  2 resistors. 


-E 


91  91 

m _ m _ 

9R4  9R32  (Ra  R4  + R32)2 


Substituting  the  partial  derivative  equivalents  in  Eqs. 
Eq.  (4)  yields 


(10) 

(ID 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(9),  (13),  (17),  and  (18)  into 
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AIm*  AKX+R4  + R32 


(RA  + R4  + R32) 


2 ARA 


<«A  + R4  + R32^ 


'RA  + R4  + R32> 


Recall  that  ARa  0.  Then,  the  relative  error  61  is'found  by  dividing  the  absolute 
error  by  the  true  value 


Therefore,  dividing  and  simplifying: 


AIm  AK  AE 

!t  ' !t  E 


RA  + R4  + R32 


RA  + R4  + R32 


The  exact  combination  of  deltas  will  not  be  known,  but  the  maximum  relative 
error  can  be  expressed  as 


AIm  ^ 

AK 

, AE 

ar4 

AH32 

Jt  ' 

Jt 

E 

RA  + R4  + R32 

RA  + R4  + R32 

To  obtain  numerical  values  for  | 6 1|  in  Eq.  (22),  thermistor  R32  was  varied  pro- 
portionally with  temperature  increasing  in  increments  of  10  degrees  Celsius. 
Values  of  all  deltas  were  taken  from  manufacturers  specifications.  It  was  calcu- 
lated from  Eq.  (1)  so  the  value  of  the  first  term  on  the  right  side  of  the  equation 
can  be  found. 

The  resistance  and  voltage  values  are  listed  on  Figure  1.  The  voltage  E is 
nominally  6.  5 volts  as  derived  from  a 6.  8 volt  ± 5 percent  zener  diode  D1  and  the 
base  — emitter  voltage  drop  of  Ql.  The  combined  error  of  D1  and  Q1  will  be  taken 
as  5 percent  for  a AE  of  0.  325  volt.  For  the  remaining  component  errors, 

AR4  0.  804  k ohms  and  AK  = 2.0  pA.  Fenwal  Electronics  thermistor  manual 
EMC-6  pages  24  and  29  characterizes  a QB41J3  standard  curve  as  listed  in 
Table  1. 
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Table  1.  Thermistor  QB41J3  Data 


True  Temperature 
<°C) 

Resistance 
(ohms  X 10'^) 

Accuracy 

. 

(±  percent  ) 

AR 

(ohms  X 10^) 

-80 

2576.  * 

14.  6 

376.  1 

-70 

1032.  * 

12.  5 

129.  0 

-60 

730.4 

10.  7 

78.  2 

-50 

389.  5 

9.  2 

35.  8 

-40 

215.  1 

7.  8 

16.  8 

-30 

123.3 

6.  6 

8.  14 

-20 

73.07 

5.4 

3.  95 

-10 

44.  76 

4.  3 

1.  93 

0 

28.  25 

3.  3 

0.  932 

10 

18.  30 

2.  2 

0.403 

20 

12.  16 

1.3 

0.  158 

Notes:  These  two  points  were  calculated  from  empirical  data  recorded  on 

21  February  1974.  The  data  is  within  +4  percent  of  the  Fenwal 
values  except  at  the  coldest  temperature  where  it  was  +6.  7 percent 
in  error. 

5jC  3$C 

See  Fenwal  Electronics  manual  EMC -6  (copyright  1974)  page  27, 
curve  No.  10  resistance  deviation  due  to  beta  tolerance  table  and 
the  explanation  below  it.  The  points  for  -80°  and  -70°  Celsius 
were  extrapolated. 

Time  Constant:  Still  air,  20  sec. 

Dissipation  Constant:  Still  air  at  25°C,  2.  5 mW/°C. 

Table  2 shows  the  contributions  of  the  terms  of  Eq.  (22)  to  | 6 1 1 . Multiplica- 
tion of  the  last  column  by  100  yields  ± percent  error  of  the  current. 

Table  3 correlates  the  relative  current  error,  degrees  Celsius  error,  and 
the  overall  percent  error  of  the  measurement.  The  degrees  Celsius  error  for 
any  temperature  was  calculated  as  the  product  of  the  following  three  terms: 
temperature  rate  of  change  with  current,  the  ideal  meter  current  It,  and  the  rela- 
tive error.  The  first  2 terms  are  from  Table  4 and  the  last  term  is  from  Table  2. 

Although  the  relative  error  in  Table  2 is  indeed  enormous,  especially  at  the 
lowest  temperatures,  it  must  be  remembered  that  this  is  the  maximum  error 
possible  with  the  given  component  tolerances.  The  major  contributor  at  the  low 
temperatures  is  the  recorder  ammeter.  A telephone  call  to  the  manufacturer 
verified  that  the  instrument's  accuracy  is  ±2  percent  of  full  scale.  The  galvanom- 
eter is  purchased  to  ±1  percent  of  full  scale  and  when  the  total  recorder  is 


Table  2.  Calculation  of  Maximum  Relative  Current  Error  Versus  Temperature 
for  the  Original  System 


Temperature 


0.  500 


0.  267 
0.  ? s8 


AE 

ar4 

~TT 

RA  + R4  + R32 

0.05 

0.  001 

0.05 

0.001 

0.05 

0.001 

0.  05 

0.  002 

0.05 

0.003 

0.  05 

0.  004 

0.  05 

0.006 

0.05 

0.006 

0.  05 

0.007 

0.  05 

0.008 

AR32 

AIm 

RA  + R4  + R32 

!t 

0.  079 
0.060 
0.043 
0.  029 
0.018 
0.  010 
0.  005 
0.  002 


Based  on  It  values  from  Table  3 (which  are  extracted  from  the  table  used  to 
process  the  recorded  data). 


Table  3.  Temperature  Error  of  the  Original  System 


Temperature 


Maximum  error 
(±  °C) 


Percent  error 
(±  percent) 


Note:  Based  on  degrees  Kelvin. 
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Table  4.  Current  Versus  Temperature  for  the  Original 
System 


Temperature 

(°C> 

Current:  It 
(pA) 

Rate 

(°C/pA) 

1 

o 

4.0 

3.0 

-60 

7.8 

2.0 

-50 

13.8 

1.  5 

-40 

23.0 

0.  85 

-30 

34.8 

0.  85 

-20 

47.0 

0.  80 

-10 

59.3 

0.80 

0 

70.0 

1.0 

10 

78.  5 

1.4 

20 

84.  6 

1.9 

considered  in  terms  of  linearity,  temperature,  function  and  ohmic  resistance,  it 
is  accurate  to  ±2  percent  FS  (full  scale). 

Determination  of  the  probability  of  any  error  in  the  envelope  defined  by  Eq. 
(22)  would  necessitate  knowledge  of  the  error  distribution  of  each  of  the  compo- 
nents. This  information  is  not  available. 

One  point  which  should  be  touched  upon  is  the  thermistor  self-heating.  Using 
the  thermistor  dissipation  constant,  listed  in  Table  1,  temperature  errors  of 
0.05°C  take  place  at  full  scale  readings,  and  significantly  smaller  errors  take 
place  at  the  temperatures  discussed  in  this  report.  Because  this  effect  is  so  much 
smaller  than  the  calculated  error,  it  has  been  ignored. 

2.  1.  2 THE  SELF  CALIBRATION 

The  manufacturer's  procedure  for  reducing  the  data  from  the  recorder  tapes 
involves  zero  and  75  pA  reference  currents.  According  to  the  recorder  manufac- 
turer, the  zero  current  indication  should  never  be  more  than  0.  25  pA  from  the 
zero  current  line  of  the  chart  paper.  Major  deviations  would  indicate  that  the 
device  is  in  need  of  repair.  The  high  current  reference  is  obtained  when  R^  and 
R32  are  switched  out  and  replaced  with  R.,  ( 82.  5 K ± 1 percent)  to  derive  a 75  pA 
current  from  the  6.  5 volt  reference  voltage.  However,  the  current  I produced  by 
applying  0.  5 volts  to  R2  is: 
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I = E/R2 

= 6.  5 volts/82.  5 k ohms 
= 78.  79  pA  . 

The  meter  has  an  accuracy  of  ±2.  7 percent  at  75  pA. 

This  calibration  is  useful  in  the  range  of  the  current  measured  but,  since  the 
majority  of  measurements  are  accomplished  at  much  lower  currents,  this  tech- 
nique is  not  valid  because  the  meter  error  K has  ever  increasing  significance  as 
current  decreases. 

2.1.3  PRECISION 

The  precision,  the  degree  of  repeatability  of  a measurement  regardless  of  its 
accuracy,  is  not  defined  in  the  flowmeter  system  literature.  However,  enough 
data  has  been  compiled  over  years  of  operational  use  to  define  the  system  precision 
as  ±0.  35  nA.  This  can  be  converted  to  degrees  Celsius  at  a specified  temperature 
by  referring  to  Table  4.  For  example,  at  -30°C  the  rate  of  change  of  temperature 
is  0.  85°C//xA.  Multiplying  by  the  precision  of  ±0.  35  ;uA  results  in  a product  of 
±0.  3°C  for  the  temperature  precision. 

2.2  The  Modified  Circuil 

2.  2.  1 ACCURACY 

The  circuit  under  discussion  is  shown  in  Figure  2.  It  is  very  similar  to  the 
standard  circuit  and  the  same  constraints  hold  for  both  circuits.  From  Ohm's  law 
it  can  be  shown  that 


Rp 


<32 


KEY 

01 

Q1 


1N4736A  68.V  + 5% 

2N1308 

4600  ohms  error  accounted  lor  by  K 
1.8k  ohms  + 1% 


^32 


16.0k  ohms  + 1% 

thermistor  UUB31J1  Ik  ohms  + 1%  et  25°C 
ammeter  accuracy  + 2%  of  full  scale 


Figure  2.  Modified  Temperature  Sensing  Circuit 

where  the  t denotes  theoretical  current  based  on  ideal  components.  If  the  compo- 
nents are  allowed  to  vary  according  to  manufacturer  specifications,  then  the  same 
expression  can  define  the  real  meter  current  for  any  set  of  components.  However, 
we  are  interested  in  the  current  as  recorded  on  the  tape.  The  measured  current 
is 




!m  = Zr  + K = E • <R4  + R32)(Ra  + Rp>  + RARP 


+ K 


(26) 


where  K is  the  meter  accuracy  as  defined  in  Figure  2 
(26) 


The  absolute  error  of  Im  is  AIm  = Im  - If.  As  a function  of  the  terms  in  Eq. 


9 1-  9I~  9Im  . aim_  . 9Im  , ^m 


^m  ^ aif  + 1#  *E  + ARA  + 3^  *RP  + TR7  *R4  *R32 

(27) 

Differentiating  Eq.  (26)  for  the  various  partials  yields 


31 

m 

9K 


= 1 . 


(28) 
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Foi  the  voltage  regulator 


9E  ' (R4  + R32)(Ea  + Rp>  + RARp 


And  for  the  resistors,  considering  the  meter  resistance. 


3Im  3 

3Ra  " E 3Ra  <R4  + R32)(Ra  + Bp)  + RARp 


E • Rp(R4  + R32  + Rp) 

[ (R4  + R32^RA  + RP^  + RARP^ 


and  the  parallel  resistance 


31  a Rp 

m tj'  3 P 

3Rp  “ * 3Rp  (R4  + R32)(Ra  + Rp)  + RARP 


E • Ra(R4  + R32) 

[ ( R 4 + R32^ra  + Rp^  + R^RpI 


91  a rp 

^ « rp  ^ _ £ 

3R4  ^ 3R4  (R4  + R32)(Ra  + RP)  + RARp 


E • Rp(RA  + Rp) 

[(R4  + R32)(Ra  + Rp)  + RARpl 


Because  of  the  similarity  in  the  positioning  of  R32  and  R4,  it  can  be  written  by 
inspection  that 


m_  E • RP(RA  + V 

*32  |(R4  + R32HRa  + Rp)  + RARpl 
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Tnis  completes  the  calculation  of  the  G partial  derivatives  for  I . The  substitu- 
tion into  Eq.  (27)  and  subsequent  division  by  1^  to  obtain  61,  the  relative  error  of 


I results  in 
m 


Alm  _AK  _AE R4  + R32  + RP 

!t  It  E (R4  + R32^Ra  + RP*  + RARP 


RA(R4  + R32)  AR 

<R4  + R32)(RA+RP)  + RARP  ' RI 


RA  + RP 

(R4  + R32)(RA  + RP»  + RARP  AI?4 


ra  + Rp 

(R4  + R32)(RA  + V + RARP  ARS2 


This  shows  that  the  absolute  errors  of  the  meter  and  the  voltage  source,  AK 
and  AE  respectively,  contribute  directly  to  the  overall  relative  error.  However, 
the  contributions  of  the  resistance  errors  are  not  so  straightforward.  Since  R^ 
was  taken  as  zero  because  it  is  covered  by  K,  the  first  resistive  term  is  zero. 
Considering  the  second  resistive  term 


RA(R4  + R32>  ARP 

(R4  + R32^RA  + RP^  + RARP  Rp 


dividing  numerator  and  denominator  by  R^iR^  + R,^)  yields 


RA  + RP 


R4  + R32 


R4  + R32 


which,  as  temperature  decreases  and 


R32  » Rp 

will  simplify  to 

1 ARp 

1 + Rp  /RA  Rp 

Considering  the  third  resistive  term 
<Ra  + Rp)  AR4 

(R4  + R32)(Ra  + Rp)  + RARP  ’ 
dividing  numerator  and  denominator  by  RA  + Rp  yields 



R4  + R32  + RaV(RA  + RP* 


(41) 


(42) 


(43) 


(44) 


which  shows  that  at  cold  temperature  where 

R32  » R4  + RaV<RA  + V ' (45) 

the  contribution  of  AR4  approaches  zero. 

The  last  term  is  similar  to  the  third.  By  inspection,  it  is  seen  that  it  will 
reduce  to 


^R32 

R4  + R32  + RARPy/(RA  + RP* 

However,  AR32  is  a function  of  temperature  and  rises  approximately  linearly. 
Even  though  the  denominator  will  become  dominated  by  R32>  Eq.  (46)  should  be 
calculated  for  a thorough  investigation  because 

R32  » R4  + RARp/<RA  + RP)  (47) 

only  below  -75°C. 
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To  calculate  the  worst  case  error,  take  the  absolute  value  of  both  sides  of 
Eq.  (37),  recalling  that  there  is  zero  contribution  from  AR^ 


AIm 

AK 

4_ 

AE 

+ -L  | 

RA(R4  + R32)  ‘ ARP 

*t 

!t 

i 

E 

RP 

+ |(ra  + Rp)  ar4| 


+ |(ra  + rp)  ar32 

where 


D - <R4  + R32^ra  + RP^  + RARP  ' 


(48) 


(49) 


Calculating  the  deltas  from  the  accuracies  of  the  components  yields 
AE  = ±0.325  volt,  ARp  =±  0,  18  k ohm  , AR4  = ± 0.  16  k ohm  and,  AK  = ±2  pA. 
The  absolute  error  of  R32  is  a function  of  temperature  (see  Table  5). 


Table  5.  Thermistor  UUB31J1  Data 


Temperature 

(°C) 

Resistance 
(ohms  X 10^) 

Accuracy 
(±  percent) 

AR 

(ohms  X 10^) 

-80 

229,  0 

10.  0 

22.9 

-70 

114.3 

7.2 

8.230 

-60 

59,  77 

5.  6 

3.347 

-50 

32,  64 

4.4 

1.436 

-40 

18.  55 

3.  4 

0.  631 

-30 

10.  92 

2.5 

0.273 

0 

CM 

1 

6.  649 

1.7 

0.  113 

-10 

4.  172 

1.  2 

0.  050 

0 

2.  691 

1.  0 

0.027 

10 

1.  779 

0.9 

0.016 

20 

1.  204 

0.8 

0.  010 

Notes:  This  table  is  based  on  data  in  Fenwal  Electronics'  Uni -Curve 
Thermistor  Manual  L-6A  (reference  number  10M-5-874-MBP) 
pages  5 and  6. 

Time  Constant:  Still  air,  10  sec  maximum 

Stirred  oil  bath,  1 sec  maximum 

Dissipation  Constant:  Still  air,  1 mW/°C. 
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Table  6 shows  the  contributions  of  the  terms  of  Eq.  (48)  to  | 6 1|  . Multiplica- 
tion of  the  last  column  by  100  yields  ± percent  error  of  current. 


Table  6.  Calculation  of  Maximum  Relative  Current  Error  Versus  Temperature 
for  the  Modified  System 


Temperature 


AK 

T 

AE 

E 

T3 

wm 

I 

0.007 

lil 

■ 

0.007 

Ksa 

0.  05 

0.  007 

mm 

0.  05 

0.007 

0.  039 

0.  05 

0.  007 

0.  031 

0.05 

0.007 

0.026 

0.05 

0.  007 

0.023 

0.  05 

0.007 

0.022 

0.05 

0.007 

0.  021 

0.  05 

0.  007 

0.  020 

0.  05 

0.007 

0.  004 
0.006 
0.  007 
0.  007 
0.  008 
0.  008 
0.  009 


0.010 
0.005 
0.  002 
0.001 
0.  001 
0.001 


0.  144 
0.  118 
0.  104 
0.095 
0.  089 
0.  088 
0.  087 
0.  087 


References 
From  Text 


Notes: 


where:  T3  = 


1 + Rp/R.  + RP/(R.  + R„) 


R4  + R32  + RaV(RA  + RP> 


R4  + R32  + RARp/(RA  + RP* 


Table  7 correlates  the  relative  current  error,  degrees  Celsius  error,  and 
overall  percent  error  based  upon  degrees  Kelvin.  The  degrees  Celsius  error  for 
a given  temperature  was  computed  as  the  product  of  the  rate  of  change,  the  ideal 
current,  and  the  relative  error.  The  first  two  terms  are  from  Table  8 and  the 
last  term  .s  from  Table  6. 


The  calculated  temperature  rise  due  to  self-heating  of  this  thermistor  under 

maximum  current  conditions,  using  the  dissipation  constant  for  still  air  listed  in 

Table  5,  is  about  0.  01°C.  This  effect  decreases  rapidly  as  the  current  decreases 

at  colder  temperatures.  Its  effect  is  so  small  that  it  will  not  be  considered. 

Resistors  R and  R . were  chosen  to  shift  the  region  of  best  accuracy  for  the 
P ** 

new  thermistor  into  the  temperature  region  where  most  of  the  data  is  gathered. 
Table  3 shows  the  least  error  occurring  between  -10°C  and  -40°C  for  the  original 
circuit,  while  Table  7 shows  that  the  modified  circuit  has  the  least  error  occurring 
between  -30°C  and  -60°C.  This  will  improve  the  accuracy  of  most  of  the  volume 


computations. 


2.2.2  PRECISION 


The  precision  of  the  modified  flow  sensor  is  the  same  as  that  of  the  original 
flow  sensor  (see  Section  2.  1.3).  In  each  case,  the  parameter  is  controlled  by  the 
meter  movement  of  the  recorder.  Precision  can  be  expressed  in  degrees  Celsius 
by  referring  to  Table  8. 

2.  2.  3 LABORATORY  TEST  AND  DATA 

A laboratory  test  of  the  modified  PR -3  flow  sensor  was  conducted  to  verify 
circuit  performance  prior  to  implementing  a balloon-borne  test.  Simulating  the 
thermistor  with  a switchable  resistance  precluded  the  errors  inherent  in  the 
thermistor  and  a temperature  controlled  chamber.  Eleven  values  of  resistance 
corresponding  to  10  degree  steps  from  -80°C  to  +20°C  were  set  with  an  impedance 
bridge  to  an  accuracy  of  better  than  0.  1 percent.  This  device  is  effectively  an 
ideal  UUB31J1  thermistor  for  the  preset  temperatures.  Table  9 shows  the  r-esults 
of  5 data  runs  with  the  ideal  thermistor  and  PR -3  serial  number  068AB. 

The  relative  error  (Table  9,  column  four)  is  well  within  the  tolerances  calcu- 
lated in  Table  6.  This  is  very  reasonable  since  all  the  data  is  gathered  using  one 
sensor  with  no  error  introduced  by  the  ideal  thermistor.  It  is  worth  noting  that 
the  precision  is  within  the  limits  mentioned  in  Section  2.  1.3.  The  next  step  was 
to  compare  the  2 thermistors  in  their  respective  circuits  during  a balloon-borne 
sampling  mission. 


¥ 


■ 


r 


Table  9.  Data  From  Laboratory  Test  of  Modified  System 


*t 

<pA) 

I 

(pA) 

AI 

Jt 

Test  Data 

Tt 

(°C) 

1 

GiA) 

_ 

2 

(pA) 

3 

(pA) 

4 

(pA) 

5 

(pA) 

— 

-80 

7.42 

7.31 

-0.015 

7.45 

7.  19 

7.45 

7.22 

7.24 

-70 

13.  89 

13.  75 

13.  7 

13.  4 

13.9 

13.  9 

13.  9 

-60 

23.  72 

23.  50 

-0.009 

23.3 

23.4 

23.4 

23.  6 

23.  7 

-50 

36.  61 

36.  50 

-0.003 

36.  6 

36.2 

36.  7 

36.  5 

36.  5 

-40 

51.00 

51.05 

0.001 

50.  9 

50.  9 

51.  1 

51.  1 

51.3 

-30 

64.  80 

65.37 

0.009 

65.  2 

65.2 

65.  5 

65.  5 

65.  5 

-20 

76.35 

77.  18 

0.011 

77.  1 

77.  1 

77.3 

77.  1 

77.3 

-10 

85.  17 

85.  94 

0.  009 

85.  8 

86.  1 

86.0 

85.8 

86.0 

0 

91.48 

92.38 

0.010 

92.  5 

92.2 

92.  5 

92.  2 

92.  5 

+ 10 

95.  85 

96.  87 

0.  Oil 

96.  8 

96.8 

97.  1 

96.  8 

96.8 

+20 

98.83 

100.  12 

0.  013 

100.  0 

100.0 

100.3 

100.  2 

100.  2 

3.  COMPARISON  OK  THE  TWO  THERMISTOR  CIRCUTS 

3.1  Experimental  Circuit 

Physical  construction  of  the  PR-3  flow  sensor  system  provides  for  monitoring 
2 thermistors.  A commutator  switches  each  thermistor  into  the  circuit,  thereby 
utilizing  the  same  voltage  regulator  and  recorder  for  each.  This  is  schematically 
drawn  in  Figure  3. 

The  similarity  between  Figure  3 and  Figures  1 and  2 is  obvious.  The  switch 
represents  the  commutator  and  the  subscript  T denotes  component?  of  the  circuit 
under  test.  If  they  are  removed,  the  original  circuit  remains.  Component  values 
are  the  same  as  in  Figure  2 for  those  with  the  T subscript  and  as  in  Figure  1 for 
the  remainder. 

3.2  Accuracy 

The  accuracy  of  the  circuit  in  Figure  3 can  be  computed  from  the  preceding 
information.  The  total  error  (±  °C)  will  be  the  sum  of  the  error  for  each  thermis- 
tor circuit,  and  will  be  computed  that  way  because  the  different  currents  produced 
in  the  2 circuits  at  the  same  temperature  preclude  the  summing  of  relative  errors. 

Since  the  same  voltage  regulator  is  used  for  both  circuits,  the  effect  of  its 
component  tolerances  should  cancel  out.  However,  since  there  is  the  difference  in 
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KEY 

R4 

60.4k  ohms +1%  j 

D1 

1N4736A,  6.8V +5% 

R4T 

16.0k  ohms +1%  “ 

Q1 

2N1308 

R32 

thermistor  QB41J3  10k  ohms  + 1%  at  25°C 

ra 

4600  ohms,  assumed  + 0%  see  K 

R32T 

thermistor  UUB31J1  1 k ohms  + 1%  at  25°C 

Rpt 

1.8k  ohms  +1% 

K 

ammeter  accuracy  + 2%  of  full  scale 

Figure  3.  Test  Circuit 


circuit  resistance  at  the  same  temperature,  it  doesn't  quite  drop  out.  The  error 
introduced  by  assuming  it  does,  is  small  and  computable.  For  example,  at  -30°C 
the  original  circuit  has  a maximum  relative  current  error  of  ±0.  153  correspond- 
ing to  ±4.  5°C,  while  the  modified  circuit  has  a maximum  error  of  ±0.  104  corre- 
sponding to  ±4.  7°C.  If  the  contribution  of  the  voltage  source  relative  accuracy, 
±0.05,  is  computed  the  relative  current  errors  are  ±1.  5°C  and  ±2.  3°C  respec- 
tively. If  the  error  in  the  voltage  supply  is  maximum,  then  the  maximum  error 
generated  by  our  assumption  is  0.  8°C.  Should  the  regulator  error  be  nonexistent, 
then  the  error  from  our  assumption  is  0°C.  The  function  is  linear  in  the  region 
between  the  examples  mentioned. 

Ammeter  accuracy  is  treated  the  same  way  since  it  is  unknown,  but  common  to 
the  2 sections  of  the  test  circuit.  In  each  case  at  -30°C,  the  meter  contributes 
approximately  30  percent  of  the  uncertainty.  Therefore,  if  the  meter  is  at  its 
maximum  error  the  error  generated  by  our  assumption  is  less  than  0.  5°C. 

Errors  for  the  remaining  components  are  as  listed  in  Table  2 and  Table  6, 
respectively.  Computation  of  the  maximum  error  °C  was  accomplished  as 
described  above.  Table  10a  is  the  error  for  the  original  circuit,  and.  Table  10b 
is  the  error  for  the  modified  circuit.  Table  10c  is  the  sum  of  the  preceding  2 tables 
and  is  the  maximum  possible  error. 
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Table  10c.  Maximum  Error  of  the 
Test  Circuit* 


Temperature 

<°C) 

Error 

(±°C) 

-70 

2.  6 

-60 

2.  5 

-50 

2.  7 

-40 

2.  2 

-30 

2.  6 

-20 

2.  7 

-10 

2.9 

0 

3.9 

10 

5.4 

20 

8.2 

Under  the  constraints  mentioned  in 
the  text. 


3.3  Kxperiim-nlal  Data 

Data  was  gathered  by  modifying  the  circuitry  in  a PR-3  Data  Recorder  so  that 
it  corresponds  to  Figure  3.  Resistors  Rprj-  and  R^-j.  were  added  to  Recorder  serial 
number  077AB,  and  R32T,  and  the  UUB31J1  thermistor  was  mounted  in  a Flow 
Sensor  next  to  the  standard  thermistor. 

Here,  the  exact  mounting  is  described.  The  flow  sensor  interior  is  radially 
divided  into  8 sections,  by  8 evenly  spaced  flow  straightening  vanes."  Mounts 
for  the  bead  thermistor  are  one  vane,  so  that  the  thermistor  leads  are  parallel  to 
the  radial  dimension  of  the  vane  and  the  bead  is  in  the  center  of  the  section  about 
0.  5 in.  (~  1. 3 cm)  from  any  surface.  The  rod  thermistor  is  secured  to  mounts 
located  along  the  circumference  of  the  cylinder  on  the  opposite  side  of  the  same 
vane.  The  mounts  are  located  so  that  the  flow  is  across  the  axis  of  the  thermistor 
which  is  also  about  a 0.  5 in.  (~  1. 3 cm)  from  any  surface.  The  distance  from 
thermistor  to  thermistor  measured  as  a cord  projected  on  the  plane  defined  by  the 
outward  facing  end  of  the  sensor  is  about  1.3  in.  (~3.2  cm).  Measuring  into  the 
sensor  from  this  plane,  the  rod  thermistor  is  back  about  1 in.  (2.  54  cm)  and  the 
bead  thermistor  is  back  about  1.7  in.  (4.3  cm).  This  equipment  was  subsequently 
used  on  3 atmospheric  sampling  missions  to  monitor  the  exhaust  port  of  a motor 
driven  sampler.  The  procedure  in  each  case  was  to  allow  the  balloon  to  rise  to  its 
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float  altitude,  then  to  activate  the  samplers  while  the  balloon  altitude  remained 
within  a nominal  envelope  to  ±1.  5 km  (5  kft).  Data  were  taken  while  air  was 
flowing  over  the  thermistors.  Temperatures  were  derived  directly  from  the 
chart  recorder  data  in  microamperes.  No  corrections  have  been  applied. 

3.3.1  FIRST  TEST  FLIGHT,  H76-7/H-78X,  13  FEBRUARY  1976 

Data  were  recorded  at  an  altitude  of  32.  2 k m ± 0.  15  k m (105.  5 k ft  ± 0.  5 k ft) 
for  280  min  (see  Table  11).  It  is  important  to  remember  that  the  data  are  com- 
mutated so  the  temperatures  listed  for  a given  time  were  recorded  about  20  sec 
apart.  Delta  (A)  is  defined  as  the  temperature  from  the  original  system  minus 
the  temperature  from  the  modified  system.  The  average  temperature  from  the 
original  system  is  -21.4°C,  and  from  the  modified  system  -22.3°C.  The  differ- 
ence of  0.  9°C  is  well  within  the  error  envelope  of  the  test  circuit. 

3.3.2  SECOND  TEST  FLIGHT,  H76-12/H-83X,  26  FEBRUARY  1976 

This  balloon  floated  at  an  altitude  of  30.  8 k m ± 0.  24  k m (101.  1 k ft  ± 0.  8 k ft) 
and  the  data  were  recorded  for  330  min  (see  Table  12).  The  average  temperature 
from  the  original  system  is  -23.3°C,  while  the  average  temperature  from  the 
modified  system  is  -20.  6 for  an  average  difference  of  -2.  6°C.  This  average  dif- 
ference is  also  within  the  computed  error  envelope. 

3.3.3  THIRD  TEST  FLIGHT,  H76-48/H-93,  25  AUGUST  1976 

Data  were  recorded  at  an  altitude  of  21.  5 k m ± 0.  03  k m (70.  8 k ft  ± 0.  1 k ft) 
for  180  min  (see  Table  13).  The  average  temperature  from  the  original  system 
is  -37.  6°C,  and  -37.  2°C  for  the  modified  system.  Once  again  the  difference  of 
-0.  4°C  is  well  within  the  allowable  envelope. 

3.1  Data  ((valuation 

Data  from  3 dual  thermistor  flights  described  in  the  preceding  paragraphs  of 
Section  3 will  be  used  to  evaluate  the  first  characteristic  mentioned  in  the  intro- 
duction; that  is,  data  exhibiting  an  oscillation  about  some  reference  temperature. 
Data  in  Tables  11,  12  and  13  definitely  have  oscillations  and  both  thermistors  track 
very  well  in  all  cases.  Tables  11  and  12  exhibit  2 1/2  cycles,  while  Table  13  has 
2 complete  cycles. 

There  is  no  common  frequency  of  oscillation.  Also,  in  each  case  the  differ- 
ence between  the  maximum  and  minimum  temperature  values  for  each  thermistor 
is  about  the  same  - that  is,  9°C,  21°C,  and  5°C  respectively  for  the  three  Rights. 
The  maximum  variation  in  temperature  due  to  the  circuit  would  have  to  be  based 
on  the  circuit  precision.  Allowing  for  ±1 /4°C  per  thermistor  circuit,  the  maximum 
difference  from  this  cause  would  be  1/2°  Celsius. 
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Table  11.  Data  From  the  First  Flight  Test 


Table  12.  Data  From  the  Second  Flight  Test 


Time 

(Min) 

Original  Part 

Modified  Part 

Current 

(pA) 

Temperature 

<°C) 

Current 

(pA) 

Temperature 

<°C) 

0 

31.  0 

-33.  0 

60.  5 

-33.3 

10 

37.0 

-28.2 

62.0 

-33.  2 

20 

41.  0 

-24.  8 

71.  0 

-24.9 

30 

47.0 

-20.  0 

77.  0 

-19.4 

40 

44.  0 

-22.  5 

76.  5 

-19.9 

50 

46.  0 

-20.8 

75.  0 

-21.3 

60 

46.0 

-20.  8 

79.0 

-17.3 

70 

49.  0 

-18.4 

78.  0 

-18.4 

80 

48.0 

-19.  2 

83.  0 

-12.  8 

90 

43.0 

-23.2 

83.  0 

-12.8 

100 

42.0 

-24.0 

77.0 

-19.4 

110 

41.0 

-24.8 

74.0 

-22.  2 

120 

43.0 

-23.2 

75.  0 

-21.3 

130 

48.  0 

-19.  2 

80.0 

-16.2 

140 

47.0 

-20.0 

77.  0 

-19.4 

150 

36.  5 

-21.  6 

73.0 

-23.  2 

160 

33.0 

-31.  5 

69.  0 

-26.  6 

170 

38.0 

-27.4 

72.0 

-24.  1 

180 

42.  0 

-24.0 

74.  5 

-22.2 

190 

45.  5 

-21.7 

76.  0 

-20.4 

200 

39.0 

-26.6 

75.0 

-21.3 

210 

35.  0 

-29.8 

69.0 

-26.  6 

220 

32.0 

-32.2 

66.  0 

-29.  1 

230 

35.  0 

-29.8 

70.  0 

-25.8 

240 

39.  0 

-26.6 

70.  5 

-24.4 

250 

40.0 

-25.7 

69.  0 

-26.  6 

260 

50.0 

-17.  6 

78.  0 

-18.4 

270 

48.0 

-19.  2 

81.0 

-15.  1 

280 

51.0 

-16.  8 

80.0 

-16.  2 

290 

58.0 

-11.0 

82.  5 

-13.4 

300 

51.  0 

-16.  8 

82.0 

-14.0 

310 

52.  5 

-15.  5 

81.  5 

-14.  5 

320 

41.0 

-24.8 

81.  0 

-15.  1 

330 

41.  0 

-24.8 

82.0 

-14.0 

A few  words  are  now  appropriate  about  the  segment  aperture  times  and  the 
thermistor  time  constants.  The  duration  of  each  segment  is  37.  5 sec.  The  time 
constants  of  the  thermistors  are  20  sec  for  the  rod  and  10  sec  for  the  bead.  Cur- 
rent recorded  for  the  2 segments  was  read  toward  the  end  of  the  first  segment  and 
the  beginning  of  the  second  segment.  This  allows  about  40  sec  separation  for  the 
temperatures  written  as  a set.  The  commutator  selects  first  the  original  (rod) 
thermistor  and  then  the  bead  thermistor.  Forty  sec  is  4 bead  thermistor  time 
constants.  In  4 time  constants,  the  bead  has  essentially  completed  any  change 
which  would  take  place  if  a new  temperature  was  impressed  upon  it  at  the  comple- 
tion of  the  segment  dedicated  to  the  rod  thermistor. 
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i.  H.iun  n\r  \ 


The  data  displayed  in  Table  14  is  from  18  flights  over  a period  of  27  months, 
and  all  the  data  is  from  PH -3  sensors  of  the  original  configuration  The  letter 
preceding  the  flight  number  identifies  the  launch  location  as  New  Mexico  (H), 
Panama  Canal  Zone  (P),  or  Alaska  (A),  Each  ambient  temperature  used  as  a con- 
trol was  obtained  from  a radiosonde  set  manufactured  by  VIZ  Manufacturing  Com- 
pany which  was  suspended  beneath  the  samplers  during  the  flight.  VIZ  Model  1091 
is  designated  AN/AMT-12  by  the  U.  S.  Army.  Specifications  published  by  the 
manufacturer  list  the  temperature  error  as  "less  than  0.45°C  average  RMS'  over 
the  range  of  +60°C  to  -90°C.  Independent  research  (Morrissey4)  cautions  that 
this  should  be  considered  nonrigorous  in  terms  of  an  absolute  figure  on  any  data 
run,  but  documents  the  "root  mean  square  deviations  in  temperature  measure- 
ments on  simultaneously  released  radiosondes  of  0.4°C." 

1.1  Data  Kwilimlion 

Data  from  Table  14  will  be  used  to  evaluate  the  second  characteristic  men- 
tioned in  the  introduction,  that  is,  the  large  difference  in  temperature  between 
similar  flights,  with  flights  at  the  same  altitude  ±10  k ft.  To  accomplish  that, 
the  average  ambient  air  temperature  on  each  flight  will  be  a control  to  which  the 
average  inlet  air  temperature  is  compared.  This  relative  temperature  as  defined 
on  Table  14  is  plotted  versus  ambient  temperature  for  each  flight  in  Figure  4.  The 
3 float  altitudes  are  differentiated  by  the  symbols  indicated  on  the  figure. 


Table  14.  Temperatures  Recorded  on  Recent  Flights 


Flight 

Number 

Date 

Nominal 
Altitude 
(k  Ft) 

Ambient 

Temperature 

(°C) 

Inlet 

Temperature 

(°C) 

sjt 

Relative 

Temperature 

<°C) 

H-29 

22  Aug  73 

70 

-54.4 

-45.7 

8.7 

H-31 

28  Aug  73 

90 

-45.  7 

-42.4 

3.3 

H-34 

6 Nov  73 

80 

-51.  7 

-44.4 

7.3 

P-138 

12  Mar  74 

80 

-59.3 

-53.  3 

6.0 

P-140 

14  Mar  74 

70 

-62.  5 

-53.  8 

8.  7 

H-37 

12  Apr  74 

80 

-53.3 

00 

1 

5.  5 

4.  Morrissey,  J.  F.  (1972)  Atmospheric  Temperature  Measurement  Using 

Balloons  and  Rockets.  AFCRL-TR-73-0380  (AD  764185),  Copyright  1972 
Instrumentation  Society  of  America. 
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The  solid  line  represents  a least  square  linear  regression  formulated  by: 

y = 11.09  - 0.36  x (50) 

where  x is  the  ambient  temperature  and  y is  the  relative  temperature,  both  in 
degrees  Celsius.  Ordinates  of  3.3°C  and  10.  5°C  were  calculated  for  abscissas  of 
-40°C  and  -60°C,  respectively,  to  locate  the  line.  The  broken  line  represents 
the  error  envelope  listed  in  Table  3.  The  empirical  data  fits  the  theoretical 
curves  very  well. 

3.  COMCI.I  SIONS 

This  analysis  began  with  a description  of  2 phenomena  which  were  to  be  vali- 
dated or  disproved  by  analyzing  the  data  system.  Errror  envelopes  for  the  cir- 
cuits involved  in  the  investigation  were  derived  and  presented  in  Tables  3,  7,  and 
10c.  Data  were  gathered  in  2 modes  to  test  the  phenomena. 

3. 1 Kirsl  Phi-noim'non 

Temperature  oscillations  during  any  1 flight  were  examined  with  the  dual 
thermistor  flights.  In  each  flight,  each  thermistor  exhibited  similar  temperature 
excursions  in  frequency  and  magnitude.  The  resistance  of  the  rod  thermistor  in 
the  temperature  range  under  discussion  ( — 30°  Celsius)  is  approximately  123  k 
ohms  while  the  bead  thermistor  is  about  Ilk  ohms.  At  these  resistances,  no 
interference  was  experienced  from  radiated  radio  frequency  energy  in  laboratory 
tests.  Therefore,  excursions  in  the  temperature  data  are  most  probably  indicating 
actual  changes  in  the  sampler  inlet  air  temperature.  Should  the  reader  be  tempted 
to  compare  the  data  from  the  two  thermistors  tabulated  in  sets,  be  cautioned 
to  remember  that  the  readings  are  actually  separated  by  a 37.  5 sec  commutator 
segment.  Considering  the  time  constants  discussed  in  Section  3.4,  it  is  reason- 
able to  expect  temperature  differences  in  the  data  "sets." 

3.2  N't'onH  1'hrnomt‘iioit 

I.arge  differences  in  average  temperatures  among  "similar"  flights  were  exam- 
ined by  use  of  existing  flight  records.  Figure  4 summarizes  the  data  from  18 
flights  and  compares  it  to  the  appropriate  error  envelope.  One  of  the  18  points 
falls  out  of  the  error  envelope  by  about  1°C  Considering  the  assumptions  made 
about  the  voltage  regulator,  ammeter,  and  radiosonde,  this  is  excellent.  It  says 

*Performed  with  Hewlett-Packard  9810A  calculator  with  the  Statistics  Block. 
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tors  were  erroneously  purchased  at  a 10  percent  resistance  tolerance  rather  than 
the  1 percent  resistance  tolerance  on  which  the  calculations  are  based.  Using 
these  10  percent  tolerance  thermistors  would  have  spread  the  data  over  a larger 
temperature  range. ) 


3.3  (ienerul  Comments 

As  mentioned  in  the  Introduction,  the  author  is  not  going  to  pursue  the  mech- 
anisms which  could  cause  temperature  changes,  but  he  will  mention  several  ideas 
which  should  be  explored.  It  is  the  author's  opinion  that  the  observed  temperature 
changes  are  occurring  in  the  air  mass  surrounding  the  balloon  system  and  are 
intimately  linked  to  the  heat  balance  of  the  entire  mass.  Energy  absorbed  by  the 
system  through  direct  sunlight  and  albedo  are  significant  factors.  However, 
waste  energy  from  the  motors  and  batteries  could  also  be  significant,  especially 
if  there  are  eddy  currents  involving  the  sampler  exhaust.  Determining  the  exact 
cause  of  the  reported  phenomenon  is  beyond  the  scope  of  this  report. 

As  a parenthetical  remark,  it  is  worth  noting  that  the  temperature  of  the  inlet 
air  is  always  above  the  ambient  temperature  reported  by  the  radiosonde.  Linear 
regressions  performed  on  the  data  in  Figure  4 in  sets  based  on  altitude  yield  a 
family  of  curves  the  slopes  of  which  increase  in  magnitude  as  the  altitude 
increases.  This  hints  that  a major  factor  in  energy  dissipation  is  convection  and 
that  as  the  air  density  drops  there  is  less  cooling.  Presently,  not  enough  data  at 
one  altitude  and  latitude  are  compiled  to  pursue  this  theory  and  it  is  beyond  the 
scope  of  this  report. 
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